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ABSTRACT 

We propose a method to improve lattice operators composed of Wilson fermions 
which allows the removal of all corrections of 0(a), including those proportional 
to the quark mass. It requires off-shell improvement of quark fields and compos- 
ite operators, which is achieved by studying the behaviour of quark and gluon 
correlation functions at large momenta. 
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1 Introduction 



In present lattice computations an important source of errors is the finiteness of the lattice 
spacing, a. A systematic method for reducing discretization errors order by order in a 
was proposed by Symanzik [||] and developed by Liischer and Weisz p[. It consists of 
modifying the action and operators by "irrelevant" terms, chosen in such a way that the 
convergence to the continuum is accelerated. In the first implementations of this procedure, 
the improvement coefficients were computed in perturbation theory, leaving errors of 0{agQ) 
in physical quantities It turns out, however, that, for lattice spacings used in present 
simulations, one-loop perturbation theory is insufficiently precise for some quantities 0, 
even when it is "tadpole improved" [^]. In refs. [|^ |^ a method for determining the 
improvement coefficients beyond perturbation theory, based on the enforcement of Ward- 
Takahashi identities (WTI's) has been proposed and implemented, thus achieving full 0{a?) 
improvement. This method allows one to determine the improved action and vector current 
for massive quarks (see also ref. p), but can only be applied to other operators, including 
the axial current, in the chiral limit. The reason for this limitation and a strategy to 
work outside the chiral limit valid in the case of on-shell matrix elements was discussed 
in ref. Vector and axial WTI's, possibly in conjunction with the use of unequal quark 
masses, have been employed to determine the coefficient of the term proportional to the mass 
in the normalization constants of certain quark bilinears |Tl] , |T2[. The question of O(amq) 



corrections is particularly significant for applications involving heavy quarks, where these 
terms may be large. Important examples in the phenomenology of S-mesons include the 
calculations of the leptonic decay constant /s, the form factors of semileptonic B decays, 
the S-parameters and the amplitudes of the radiative decay B K*'y. 

In this paper we suggest an alternative method to extend operator improvement beyond 
the chiral limit, based on the consideration of the high momentum behaviour of quark 
correlators. The basic ingredient is the matching of lattice and continuum Green functions 
at large Euclidean momenta (or equivalently at short distances), where perturbation theory 
applies. An important property, crucial for our purposes, is that renormalized continuum 
correlators do not contain contact terms beyond those found in perturbation theory. This is 
a general consequence of the Callan-Symanzik equation, when applied at large momenta in 
an asymptotically free field theory. We stress that we aim only to remove errors proportional 
to a and do not attempt to work at arbitrary values of arriq, as in the program proposed in 
ref. |T^. In particular, after implementation of our proposal there will still remain errors 



proportional to a^m^. 

Our method allows for the evaluation of the renormalization constants of divergent operators, 
such as the scalar and pseudoscalar densities. For a discussion of recent determinations of 
Zp, Zs and their ratio see ref. In principle this approach can be applied to composite 
operators of arbitrary dimension, though at the price of a very high proliferation of mixing 
operators, which may render the actual implementation of this scheme impractical. 



^By niq we mean some choice of physical quark mass. We do not need to pick a specific definition in this 
paper. 
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We wish to mention, as a side remark, an important limitation of the general idea of ex- 
ploiting the large-momentum behaviour of correlators to renormalize composite operators 
non-perturbatively. The point is that for composite operators which can mix with lower 
dimensional ones, such as the octet part of the weak Hamiltonian, the strategy of killing at 
large momenta the "form factors" not allowed by chiral symmetry cannot be made to work. 
In fact lower dimensional operators may give rise to a large-momentum behaviour that can- 
not be distinguished from that due to the exchange of Goldstone bosons. For completeness 
we will briefly discuss this question in Appendix A. 

An alternative method [^, which does not suffer from these limitations, is to measure on the 
lattice the behaviour of correlation functions in x-space in the region a <^ \x\ <^ ^qcd^ where 
perturbation theory is reliable, and use the O.P.E. with perturbatively computed Wilson 
coefficients to fit this dependence. The parameters determined by the fitting procedure are 
the physical matrix elements of the (finite, renormalized) operators that appear in the O.P.E. 
In this way there is no need to know the expression of the renormalized operators in terms 
of the bare operators of the regularized (lattice) theory. The question here is whether the 
lattice is sufficiently fine-grained for the window a -C \x\ -C ^qcd exist. If this is the case, 
the method can be applied to a number of interesting cases, some of which are discussed in 
refs. |15| and ||16| . 



The approach presented in this paper is a natural extension, valid up to O(a^), of the method 
proposed in ref . , where the renormalization procedure is carried out in a non-perturbative 
way on correlators with external quarks and gluons. As discussed in ref. , such a procedure 



is expected to work if we can choose the virtuality of the external states, /i, at which the 
normalization is carried out, so as to satisfy Aqcd -C -C 1/a. Without improvement, the 
relative errors in matrix elements are of order arrig and oAqcd- In principle, our method is 
able to remove these errors and leave only corrections of O(a^). 

As has been said above, the condition Aqcd <^ /i is required so that one can use perturbation 
theory to relate matrix elements renormalized on the lattice to those computed in standard 
continuum schemes. We wish to stress that, in principle, we can relax this constraint by 
using a sequence of lattices of decreasing a and correspondingly decreasing physical volumes. 
The idea is to determine the normalization constants by imposing renormalization conditions 
on a lattice with a very small lattice spacing. On such a fine-grained lattice large values of 
H are allowed, but the smallness of the physical volume precludes the possibility of reliably 
computing hadronic correlation functions. It is therefore necessary to determine the corre- 
sponding normalization constants for the larger, but coarser, lattice on which the physical 
matrix elements are finally computed. This is achieved by increasing the lattice spacing 
a and matching the normalization constants at fixed /x, and successively decreasing and 
matching at fixed a. The procedure is repeated until the required value of a is reached. This 
method, which has been proposed within the context of the Schrodinger functional in ref. , 
could be applied also here. A first attempt in this direction has been made recently |T3 . 



Our approach requires fixing the gauge, a procedure usually afflicted by Gribov ambigui- 
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ties [0 0. However, in the large-momentum region, where correlators are well described by 
perturbation theory, these ambiguities vanish asymptotically. Gauge fixing does however, 
introduce a further complication. When we evaluate correlators of gauge invariant operators 
with external quark and gluon fields, we need, in general, to include mixing with gauge non- 
invariant operators. The form of the latter is restricted by BRST symmetry [^, For 
flavour non-singlet quark bilinears, these restrictions are sufficient to forbid the appearance 
of gauge non- invariant terms until O(a^). For higher dimension operators, such as the four- 
fermion operators in the effective weak Hamiltonian, gauge non-invariant operators should 
be considered even at 0{a^). The investigation of the consequences of BRST symmetry for 
the evaluation of Green functions of gauge-invariant operators between off-shell quark and 
gluon states is one of the principal goals of this paper. Our conclusions and the imphcations 
for the improvement of lattice actions and operators are presented in Sec. ^ below. 

The paper is organized as follows. In Sec. |] we explain our strategy and apply it to the im- 
provement of quark flelds. In Sec. ^ we discuss the procedure to improve quark bilinears fully 
up to O(a^). Few conclusive comments can be found in Sec. ^. The paper ends with three 
Appendices. In Appendix A we discuss in a simple example why in presence of spontaneous 
chiral symmetry breaking it is not possible to renormalize non-perturbatively operators that 
can mix with lower dimensional operators by only looking at the high momentum behaviour 
of correlators. In Appendix B we give a simple argument to prove that in off-shell ampli- 
tudes BRST-symmetry does not forbid the mixing of gauge-invariant operators with gauge 
non-invariant ones, if the latter vanish by the equations of motion. In Appendix C we apply 
our off-shell improvement procedure to the quark propagator calculated at one-loop order in 
perturbation theory. 



2 Strategy and implementation on bilinears 

Following ref . |^ , we recall the basic steps of the Symanzik improvement program for lattice 
QCD up to O(a^). We will refer to the standard Wilson formulation of the theory. The 
improvement of the spectrum can be accomplished by adding the SW-Clover operator to the 
lattice QCD action with a coefficient determined using suitable WTI's identities for the axial 
current. This coefficient is a function of g^, and has been determined non-perturbatively for 
several values of Qq using numerical simulations [§. The improvement of the action is not 
sufficient to remove 0(a) terms from correlation functions ^. To achieve this result one 
must also improve the relevant operators, by adding to them terms of three kinds 

(a) Operators which vanish by the equations of motion, which can be either gauge invariant 
or non- invariant. They give rise only to contact terms, and thus do not affect the long- 

■^For recent progress on this issue see refs. 

■*This classification was introduced in the study of rcnormahzation of gauge invariant operators, which 
requires, in general, operators of all three types ||2^ The classification applies as well for improvement, 
since the allowed operators are restricted by the same type of symmetries relevant for mixing. 
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distance properties of correlators. 



(b) Gauge non-invariant, but BRST invariant, operators, which also do not contribute to 
physical amplitudes. 

(c) Gauge-invariant operators which do not vanish by the equations of motion. These are 
needed to cancel discretization errors proportional to a in physical amplitudes. 

Within the approach of ref. 0, by working in the chiral limit, one is able to improve operators 
by using only physical matrix elements. Improvement of composite operators in on-shell 
matrix elements was extended beyond the chiral limit in ref. |jlO| (see also refs. |1TI| , |12|)- In 



both cases only terms of type (c) are required. Here we need to consider all three types of 
terms. 

The idea developed in this paper is to construct finite operators with well defined chiral 
transformation properties by matching lattice and continuum correlators in the region of 
large Euclidean momenta where perturbation theory applies. This is most easily done by 
requiring the vanishing of chirality-violating form factors in correlation functions in which 
the operator is inserted together with external quark and/or gluon (and possibly ghost) legs 
with large virtualities. This procedure is justified by the following two observations. First, 
at large momenta renormalized perturbation theory becomes chirally invariant (note that: 
i) explicit chiral symmetry breaking effects induced by the regularization can be reabsorbed 
by imposing the validity of the WTI's of chiral symmetry and ii) violations from the non- 
vanishing of quark masses disappear at large momenta). Secondly, contributions due to the 
spontaneous breaking of chiral symmetry, which are absent in perturbation theory, die off 
at large momenta. So the magnitude of both effects (which come from both perturbative 
and non-perturbative violations of chirality) decreases as we go deeper into the Euclidean 
region. As discussed in Appendix A, the argument cannot be applied to determine the 
mixing coefficients of lower dimensional operators. 

In the following we consider correlation functions of elementary fields and multiplicatively 
renormalizable operators. They have, in general, non-vanishing anomalous dimensions. To 
discuss the improvement of such operators in a systematic way, one must first multiply them 
by appropriate renormalization constants so that they have a finite continuum limit. One 
can then speed up the approach to the continuum limit by adding further operators of 0(a), 
and by allowing the normalization constants of fields and operators to depend on aniq. Thus, 
to the order at which we work, all the Z-factors we introduce will depend linearly on amq. 
We do not, however, exhibit this dependence explicitly - it is fixed automatically by the 
overall normalization conditions that will be imposed. 



2.1 Lattice BRST symmetry 

In order to classify operators of type (b), i.e. gauge non-invariant operators which are BRST 
invariant, we need to know the form of the lattice BRST transformation. This is a simple 
generalization of the continuum transformation, as explained in ref. pi|. 
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The lattice Landau gauge can be enforced by requiring that 



G(f/)=ETr(f/„,, + f/i,) (1) 



n,fi 



be locally maximized along a gauge orbit. This is equivalent to imposing, at each site n, and 
for each color a, 



^ m) ^ ^«<^') 



I 



u;=0 ^ 



Tr \t\\J^^, - Ul;) - f^iU^.,,, - Ul_,,,)] , (2) 



where u;^ parameterizes the gauge transformation at site n and Qn = exp(ia;^t"). The 
condition (0) is the discretized version of (9^Aj^(x) = 0. We can implement the lattice 
Landau gauge in the functional integral in the standard way by adding to the action gauge 
fixing and ghost terms of the form 



^9f = «^ E 



-aA„A„ + iXJ,XU) + 2^ c„ ^^^^ 



n' 



a;=0 



(3) 



In eq. (|^) c and c are ghost and anti-ghost fields, while A is a Lagrange multiplier, all of 
which are defined on lattice sites. We obtain the lattice Landau gauge in the limit a — 0, 
for then the integral over A^ sets = 0. 

There are, in general, many solutions to the gauge condition Eq. (0): this is the problem 
of Gribov copies As mentioned in the introduction, we assume here the absence of 

copies for the gauge configurations which are responsible for the dominant contributions to 
correlation functions at large momenta. 

The lattice action, 5, and the gauge- fixing action, iS^/, are separately invariant under the 
appropriate BRST transformations. The latter are the lattice generalization of the contin- 
uum transformations. They are constructed by letting gauge and fermion fields be gauge 
transformed with a gauge matrix 5'„(c) = exp(iec^t°), where e is a constant Grassmann 
parameter: 

e 5brst q-n = gn{c)qn - qn = iec'^fgn (4) 
e^BRSTgn = -iqneCt" (5) 
e 5brst t/n,M = (<^"f/„,M - UnJ^c^^^) . (6) 

Note that higher order terms vanish since = 0. The other fields transform in the usual 
way, i.e. 

^BRST = —-fabcCnCn , (^BRST C° = A^ , (5brsT -^^ = . (7) 

The lattice BRST transformation is nilpotent, like its continuum counterpart. 

The full action is also invariant under the anti-ghost shift symmetry, ^ + const. This 
symmetry is also very useful in restricting the form of improved operators. 
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2.2 Improvement of the action 



The Wilson fermion action is improved by the addition of the SW operator, 

<S = <S gauge + '^Wilson ~^ J '^'^^ ^SW (8) 

= ^gaugc + J d'^x q0 + mo)q , (9) 

where 

Osw = -jcsw ^qcTf^yF^^q, (10) 

with a^j^u = (i/2)[7^, 7i,]. In Eq. (y), + mo is a shorthand for the entire lattice fermion 
operator. Here and in the following we use continuum notation to refer to lattice quantities. 
For physical quantities to approach their continuum values with deviations of only O(a^), 
we need, in addition to choosing the correct value of csw, to adjust the bare gauge coupling 
go and the bare mass tuq appropriately, in a way which depends on the renormalized quark 
mass. This is discussed in detail in ref. [0, but will not concern us here, since we have in 
mind fixing a{go) and mo using hadronic quantities, such as m^ and m^^. 

There are two additional gauge invariant operators of dimension 5 which could, in principle, 
be added to the action. The first one is 

S[ = arriQc'i j d^xqi^ + mo)q , (11) 

but this can be immediately eliminated by rescaling the quark fields. The second one can 
be written in the form 

S'^ = 2ac'2 J d^xq(P + mofq (12) 

and thus vanishes by the equations of motion. Although it can be neglected for the compu- 
tation of on-shell quantities, it gives rise to contact terms in correlators. Its contributions 
can, however, be reabsorbed by a suitable redefinition of the 0(a) part of fundamental fields 
and operators, as discussed below. Consequently we do not include it in the action. 

Finally, we must consider the possibility that gauge non-invariant operators may need to be 
added to the action in order to improve quark and gluon correlation functions. Any such 
operator should be invariant under both the lattice BRST symmetry and the anti-ghost shift 
symmetry, and also under global color transformations, lattice rotations and translations. In 
addition, it must have zero ghost number. Allowed terms thus have the form (again using 
continuum notation for lattice quantities) 

Sgni = j d^x 5brst [c'^^m^^] ' = X'^id^c, c. A, q, q, A) , (13) 

where the vector is a color octet with zero ghost number. The possibility of lowest 
dimension is X" = A" with some choice for the lattice gauge field, ^4°. This leads, however. 
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to a term proportional to Sgf{a = 0) + O(a^), which can be absorbed by rescahng A and a. 
Improvement terms of 0(a) would result from choices of an of dimension 2. There exist, 
however, no such operators, and therefore no possible gauge non-invariant improvement 
terms at 0(a). At next order, by contrast, there exist a number of possibilities, e.g. X° = 
qi^.f'q or X^J = fabcC^d^.c''. 

2.3 Improvement of quark fields 

To improve quark fields to O(a^) we must add all possible operators of dimension 5/2 which 
are allowed by the unbroken lattice symmetries. In particular, these operators must have the 
same properties as the quark fields under global gauge transformations, rotations and flavour 
transformations, and must satisfy the same BRST identities. The bare operators having these 
properties are pq, rrioq and ^q. The appearance of ^q is, at first sight, surprising since it 
transforms differently from the quark field q under local gauge transformations [and thus 
under BRST transformations, Eq. (^)]. Nevertheless, as we now explain, the non-linearity 
of BRST symmetry does not exclude such an operator. 

The constraints which follow from BRST symmetry can be obtained by requiring the im- 
proved operators to satisfy (up to 0(a^)) the same identities that are satisfied by the con- 
tinuum operators. The continuum identities take the form 

= (<5brst[Oi(xi)02(x2)...]) 

= (fe[Oi(xi)]02(x2) . . .) ± (Oi(a;i)5B[02(a;2)] ...) + ..., (14) 

where Oi are local composite operators located at arbitrary positions. The choice of sign on 
the second line depends on whether Oi is bosonic or fermionic. From now on to lighten the 
notation we will use the abbreviation (5b for ^brst- Note that for the identities to be non- 
trivial the product of operators being varied must have ghost number —1. The conditions 
on the lattice operators become 

([5i3i]02 . . .) ± m5^2] ...) + ... = 0(a2) , (15) 

where Oi and [^bOj] are, respectively, the improved lattice versions of the continuum operator 
and of its continuum BRST variation. For brevity, we have subsumed the site label into the 
definition of the operator. The expectation value is to be taken with respect to the improved 
lattice action. 

Since the improvement of the action does not break the lattice BRST symmetry, one can 
derive analogous BRST identities directly on the lattice 

(5b[Oi]02 . . .) ± {dM02] ...) + ... = , (16) 

where the variations are under the lattice transformations (§)-(0). These identities are of 
the required form ([ISD if 

[5^^]=5^[^i\ + 0{^). (17) 
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In other words, the improved lattice version of the continuum BRST variation on the l.h.s. 
must equal (up to 0(a^)) the lattice BRST variation of the improved operator on the r.h.s. 

We now discuss the consequences of the condition (|T^. Since lattice and continuum BRST 
symmetries take the same form, ([T7|) is automatically satisfied up to 0(a). Eq. (0) is thus 



a constraint on the form of improvement terms. 

Consider first a gauge-invariant operator composed of quark and gluon fields. For this the 
continuum BRST variation vanishes and the condition (|l^) simply becomes (5b [Oj] = O(a^). 
In other words, the terms of 0(a) added to improve Oj should themselves be invariant un- 
der the lattice BRST transformation, with the exception of operators which vanish by the 
equations of motion. These do not need to be BRST invariant because they contribute 
only when the insertion points of two operators happen to coincide, in which case the con- 
straint (|T7p must be applied to the resulting composite operator. We do not discuss this 
issue further since, as explained below, BRST non-invariant operators are not needed for 
O(a^) improvement. 

For the case of a quark field the continuum BRST variation is non-vanishing and the condi- 
tion ([T7| ) becomes 

ie[^q] = 5B[q]+0{a^). (18) 

As noted above, symmetries other than BRST require that the improved quark field has the 
form 

Here Zg contains an implicit dependence on the quark mass, and we have chosen to group 
^ with niQ since this combination gives only contact terms in correlation functions. From 
Eq. (H) the lattice BRST variation of the improved quark field is 

5b [g] = {c*^ [l + ac'gC^ + mo)] q + acNci ^[c^q]} ■ (20) 

Thus, comparing (|^) with ( plSf ) and using (p!9D , we learn that 

= c^q + aZ;'/^c^Gi[n ■ (21) 

Naively, we might have expected the second term to be absent, since the ghost field itself 
does not require improvement at 0(a). But in fact we have no a priori knowledge of how to 
improve the composite operator c"g, and there is no inconsistency with having the additional 
term proportional to cngi- We conclude that BRST symmetry does not forbid non gauge 
invariant improvement terms for the quark field. For completeness we note that the improved 
antiquark field takes the form 



q = Zg q\l + acg(-^ + mo) - acNci' 



:22i 



Having determined the general form of the improved quark and antiquark fields, we now 
explain our method for determining non-perturbatively the improvement coefficients c'g and 
Cngi and the normalization factor, Zg. 
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To the order in a at which we are working, Zq has a hnear dependence on the quark mass, 
while c'g and cngi are mass independent. The improvement condition is that the lattice and 
continuum renormalized quark propagators, S and S respectively, differ only at O(a^), i.e. 

S{p) = S{p) + 0{a') . (23) 

Here, the renormalized lattice quark propagator in momentum space is defined by 

S{p) = J d^xe-'P-''S{x) , (24) 

where 

S{x) ^ {q{x)m) (25) 
1 + ac'q{^ + mo) + acNGi ^ C[iO)q{x) 1 + acq{-ip + mo) - ocngi ^ 



In the following, we denote the bare lattice quark propagator by 5*^ {Si,{x) = {(l{x)q{0))) . 
The continuum renormalized quark propagator has the decomposition 

Sip)=tE,{p')j^+j:2{p') (26) 
with some choice of normalization condition, to be specified below. The large p"^ behaviour 



of Si and E2 is determined by the renormalization group equation |2g]. In particular, 

^i(p^) — ^ computable logarithmic corrections, while, in the chiral limit, 

^2(p^) {qq)/p'^, again up to logarithms, where (qq) is the usual order parameter for 

spontaneous chiral symmetry breaking. Away from the chiral limit, S2 still vanishes at large 
p^, but now decreases as rriq/p'^. As already indicated in Eq. (^3]), our improvement condition 
will be that the lattice propagator has the same asymptotic behaviour as the continuum one 
up to O(a^). 

The effect of the improvement terms can be seen by deriving the expression for the improved 
propagator in momentum space. One gets 

S{p) = Z;'[SLip) + 2ac'q + ac^Gl{^^^,SL{p)} + 0{a^)) (27) 
= (SLip) + 2ac'q) + 2acNGi \2t ^2{p') - p'Si(p2)j + 0{a^) . (28) 



In the second line we have used Eq. (|26| ) to parametrize the form of the propagator. The term 
proportional to c'q thus adds a constant to the propagator, corresponding to a delta-function 
contact term in position space. The contribution from the term proportional to cngi is more 
complicated. In the chiral limit, and for large p"^, we can ignore S2, and the cngi term is 
proportional to p^Si(p^). As already noted, this combination varies only logarithmically 
with p^. Thus, in the chiral limit, the effect of the two improvement terms is similar. In 



practice this leads to problems in separately determining the two coefficients c'q and cngi HI. 
Away from the chiral limit, the cngi term also contributes a mass-dependent renormalization 
of the coefficient of j/. 
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From this analysis we thus expect that the unimproved lattice propagator does not decrease 
at large momenta, but instead contains a constant and a slowly varying term which must 
be cancelled by appropriate choices of c'q and cngi- and cngi are then determined by 
requiring that the improved propagator satisfies, for a range of asymptotic values of the 
condition 

TiSip) = , (29) 

where the trace is over spin and colour indices. As Zg is just an overall factor in this 
condition, we do not need to know it in advance in order to determine c'g and cngi- The 
implementation of Eq. (p9D requires that the asymptotic behaviour of the correlation function 
can be extracted in the region where <^ 1/a^ to avoid distortions due to O(a^) errors. The 
same caveat applies to the improvement of bilinear operators discussed below. Exploratory 



numerical studies of the condition (E9f) can be found in 14 



Assuming that we have fixed c'g and cnqi in the way we have described, we can determine 
Zg, by imposing, for instance, the condition 




1 . (30) 



In view of the above discussion it may be preferable to determine Zg using the forward 



matrix element of the conserved lattice current on quark states ^7j. We note that, while 
different definitions of the field renormalization Zg are possible, the coefficients c'g and cngi 
are completely fixed by the form of the action, and are functions only of Qq. This is why 
they are completely determined by the asymptotic behaviour of S2. 

We can now see why c'2 can be set to zero in the action. The reason is that the term 
proportional to c'2, eq. (|12|), can be removed (up to 0(a^)) by performing the change of 
variables 

q ^ 1 + ac2(^ + mo) q , q ^q 1 + ac2(-^ + mo) . (31) 

The expressions for the improved quark fields in terms of the new variables retain the same 
form, Eqs. (|l^) and (|2^), but with a shifted improvement coefficient c'g — * c'g — d^- The 



Jacobian for the change of variables shifts the g\ coefficient in front of the gluon action, 
by a function of C2 [whose expansion begins at ©(cg')]. Since we are determining dg and a(5'o) 
non-perturbatively, these shifts do not concern us. Similar shifts occur in the improvement 
coefficients, c'o, of other quark bilinears which we will introduce below. 

We have implemented the off-shell improvement procedure described in this sector in per- 
turbation theory, using the one-loop results of Ref. Details are given in Appendix C. 
We find that a one-loop calculation can only separately determine the tree-level values of 
dg and cj^gi (the results are dg = —1/4 and cngi = 0). At one-loop accuracy, only the 
linear combination c'g + cngi can be determined. Thus, at present, the separate one-loop 
corrections to dg and cjsgi are not known. 
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3 Improvement of bilinear quark operators 



In this section we discuss the improvement of local gauge-invariant bilinears. We consider 
only flavour non-singlet bilinears, but we drop flavour indices because they are not necessary 
for our discussion. We introduce the following notation: 

• The amputated vertex function is deflned as in the continuum | 

(p|Orb')an.p = S{p)-'Gr{p,p')S{pr' , (32) 

where 

Grip,p') = I d'xd'ye-'^-^+"P' y{q{x)Orm{y)) ■ (33) 

with Oy{x) = q{x)rq{x) and F a Dirac matrix. Note that Or is deflned in terms of bare 
quark and antiquark flelds. We do not construct the operator using improved quark 
flelds q, because these contain the gauge non-invariant improvement term proportional 
to cngi which in the end, as we shall see, does not appear in bilinears. 

• Improvement of bilinears requires consideration of operators of the form 

Er = q[T(^ + mo) + {-'^ + mo)T]q. (34) 

These operators vanish by the equations of motion and therefore only contribute con- 
tact terms when inserted in correlators. 

It turns out that the improvement of bilinears up to O(a^) does not require the addition of 
gauge non-invariant operators. This is because there are no such operators with the same 
quantum numbers as the bilinears, having dimension 4. In fact, as mentioned above, allowed 
gauge non-invariant operators are of two types: those vanishing by the equations of motion, 
which need not be BRST invariant, and BRST invariant operators, which need not vanish 
by the equations of motion. The former appear flrst at dimension 5, an example being 

f [r7^tM^(:^ + mo) + (-^ + mo)t"A;:7^r] q. (35) 

The lowest dimension BRST invariant operator has, instead, dimension 6: 

<5b (g[r, i.rqd.c'') = q[r, 7jt'^g [d^X'^ + faM{d,c')c''] . (36) 

Our restriction to flavour non-singlet operators is crucial here. Improvement of flavour singlet 
bilinears requires, in general, gauge non-invariant operators of dimension 4. For example, 
for r = 1, one requires operators of the form appearing in the gauge-flxing action, Sgf. 



In ref. |17 the amputated vertex function was denoted Aor(p,p')- 
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3.1 Pseudoscalar density 



The renormalized improved pseudoscalar density has the general form 

P{x) = ZpP{x), (37) 
Pix) = [q{xh,qix) + ac'pE^,{x)] . (38) 

The form of the terms in the amputated vertex which are cancelled by an appropriate choice 
of c'p can be found by considering the contribution proportional to c'p itself. This can be 
obtained by calculating the contribution of E^^ to the vertex function Gp (see Eq. (|3^) ), 
with the result Q 

ac'p[Sip)^, + ^,S{p')]+0{a^). (39) 

After amputation, one finds 

(p|ac'pE^,|p%^p = ac'p [i,S{p')-' + S{pr'-f,] + 0{a') . (40) 

For large p^, p'"^ and {p — p')"^ we expect the continuum vertex function to approach its 
perturbative form, and thus that the only surviving form factor is that proportional to 75. 
Indeed, because of chiral symmetry, the form factor proportional to j^j^ vanishes in contin- 
uum perturbation theory when nig = 0. If improvement were not implemented however, a 
contribution of 0(a) proportional to a{p — p')fj,'jfj,'y5 would survive, as can be seen by sub- 
stituting S{p)^^ into Eq. (^). This may be cancelled by tuning c'p. One can thus 
determine c'p by setting p ^ p' and imposing the condition 

lim„ Tr(7^75(p|P|p')amp) = 0. (41) 

An incorrect choice of c'p would result in an 0(a) contribution to this trace growing like p^. 
Eq. (^) is equivalent to determining c'p from 

/ 1. Tr (7^75(p|g75g|p%i„p) 
— acp = hm —— — — — j — : -, (42) 

IpUp'Mp-P'Hoo Tr (7/^75 (P I ^75 Ip ) amp) 

which shows that the overall normalization of the operator need not be fixed before deter- 
mining c'p. 

The normalization constant Zp is logarithmically divergent as a — >^ 0, and should be fixed 
by a renormalization condition. A possible choice is 



^PJ^Tr (75(p|P|p)amp) 



1 . (43) 



We note that it is necessary to determine c'p before Zp. The point is that Eq. (^) contains 
a term of 0(a) which is proportional to 75, coming from the "mass term" in the inverse 
propagator, i.e. from the part of 5*"^ proportional to the identity. Tuning c'p thus changes 
the 75 form factor by terms of 0(a), and so changes the resulting value of Zp. Note also 
that an error in the determination of c'p leads to an error in the 0(a) part of Zp that is 
proportional to arriq. Similar comments apply to other bilinears. 

^We stress that Gp contains the insertion of the operator P and not P. 
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3.2 Vector Current 



The improved vector current has the general form 

%{x) = ZvV,{x), (44) 
V^{x) = q{x)-f^q{x) + acy^idy \q{x)a f,yq{x)] + acvE^^{x) . (45) 

V 

In this case, improvement requires the determination of three constants, one more than for 
the pseudoscalar density. The new type of constant, cy, is to be chosen to remove 0(a) 
contributions to physical amplitudes in the chiral limit. It was introduced and discussed in 
ref. and determined in refs. [|12], The contribution proportional to cy is needed to 
cancel contact terms of 0(a). Finally, the part of Zy proportional to aruq (usually called hy) 
should be chosen so that the vector current is normalized as in the continuum up to O(a^). 
The normalization Zy is not arbitrary, unlike that of the pseudoscalar density, but should 
be fixed so that the vector current satisfies the usual current algebra relations. Zy is finite 
in the continuum limit. 

To determine the improvement constants cy and dy we use the fact that, for the vector 
current, the only continuum form factor which survives at large momenta is that proportional 
to 7^. The contribution of c'y to the amputated vertex, (p| V^|p')amp, is 

ac'y [i^.S{p')-^ + 5(p)-Sj + 0(a2) , (46) 

which for large \p'\ and |p — p'| is proportional to 

acV [7/. /+ vhA + 0(a2) = ady [{p + p')^ + ia^,{p - p'),] + 0{a^) . (47) 

The contribution of cy, on the other hand, is asymptotically proportional to a^y{p — p')u- 
In the unimproved current, therefore, contributions of 0(a) proportional to the identity and 
to cr^jy are not suppressed, and must be cancelled by tuning cy and c'y. We can determine 
c'y separately by setting p = p', since in this case the cy term does not contribute, being 
proportional to a total divergence. To isolate the unwanted form factor we project onto the 
identity matrix, and impose the asymptotic condition 

lim Tr((p|V;b)a^p) = 0. (48) 

>oo 

Having determined c'y, we then fix cy by requiring that the tensor form factor is absent 
asymptotically, i.e. 

TtMpK\p'),^,^) = 0. (49) 

Here and in the following, there is no implied summation over repeated indices unless ex- 
plicitly indicated. Since the tuning of c'y affects both scalar and tensor form factors, it is 
necessary that the conditions ( ^8]) and (^91) be applied in the specified order. 
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Finally, the normalization constant Zy can be obtained by imposing the WTI 

E(pl^M^Mb')amp = S{v)-^ - SipT' , (50) 

where 9^ is the symmetric (improved) lattice derivative 

^^,f{x) = ^ [/{x + afi) - /(x-a/2)] . (51) 
, this gives the same result as the condition 



= 1 . (52) 

Note that, if we express all improved quark fields in terms of bare quark fields, using Eqs. ([T9|) 
and (p2D, both sides of Eq. (|^) will be proportional to Zg. Thus we do not to need to know 
Zg in order to determine Zy using the WTI. The factors of Zg simply ensure that each term 
in Eq. remains finite in the continuum limit. 

An alternative way of determining Zy is to require that the charge of a particular hadron 
equals its physical value. In practice, this may be preferable to the methods described above, 
which are based on quark correlators, since determinations using hadronic states typically 
have smaller errors. 



As discussed in ref. \TJ 



amp J 



3.3 Axial Current 

Improvement of the axial current is accomplished following similar steps to those needed for 
the vector current. The improved form of the axial current is 

A^{x) = ZaA^{x), (53) 
^fii^) = g(a;)7^75g(x) + acA 9^ [g(x)75g(x)] + acA-E^^ 75 (x) . (54) 

In this case, the contribution of -^7^75 to the amputated vertex, {p\A^\p')i^^p, is 

ac'A [l,l,Sip')-' + Sip)-\^,] + 0(a2) . (55) 

The only form factor which must survive asymptotically in the continuum is that proportional 
to 7^75. We therefore first determine c'a from the condition 



lim Tr (cr^^75(p|^Mb)amp) = (56) 

|p|— >oo 



and then fix ca by imposing 



lim„ Tr(75(p|A^b')amp) = 0. (57) 
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Finally, the normalization constant Za can be determined by enforcing the WTI 

For p = p' this identity fixes the value of m^, since the first term on the l.h.s. vanishes. Here, 
fhq is a possible definition of the renormalized quark mass, which is automatically improved 
to O(a^), being determined from quantities that are already improved at this order. Note 
that the combination rhqZp is determined unambiguously, but that, as in the continuum, 
the value of rhq depends on the normalization condition chosen for Zp. Having found m^, 
Za can then be obtained by considering p ^ p', and tracing Eq. ( [58| ) with 7^75. As for Zy, 
the determinations of and Za do not require knowledge of Zq. 

3.4 Scalar and tensor bilinears 

Finally, we briefly sketch the analysis for the scalar and tensor operators. The improved 
scalar density is 

S{x) = ZsS{x), (59) 
S{x) = q{x)q{x) + ac'sEsix) . (60) 

We can fix c's from the condition 

lim Tr(7^(p|5b)amp) = (61) 

|p|— >CXD 



and then determine Zs by imposing 

1 

T2 



^S'TT^Tr ((p|S'|p)amp) 



1 . (62) 



The improved tensor operator has the form 

T^y{x) = ZrTf^^ix) , (63) 

T^u{x) = iq{x)a^yq{x) + acT(<9^K - dyVf,) + ac'rET^^ix) , (64) 

where in V^, defined in Eq. (^), one does not need to include the 0{a) terms proportional 
to Cy and c'y- We can obtain c't by requiring 

lim Tr (7^75 (p|T^i.b) amp) = , (65) 

where all Lorentz indices are different, and then determine ct by enforcing the condition 

lim Tr (7A.(p|7].i.b')amp) = . (66) 

\p\ ,\p \ ,\p—p \-^QO 



Finally, we can fix Zt by imposing 

1 
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Zt — Tt {(r^p{p\T^p\p) ^rnp) 



1 . (67) 
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4 Conclusions 



We have presented a method for determining the off-shell improvement coefficients of gauge 
invariant bilinears in the massive case, using quark correlation functions. We stress that 
our ultimate interest is in the calculation of the coefficients necessary for the improvement 
of on-shell quantities, but that we might wish to use off-shell correlators as an intermediate 
step. A detailed numerical study will be necessary to determine whether the method is 
practical. The pilot study of ref. Jl^ suggests that it will be difficult to separate c'q from 
Cj^Gii in which case the method works only in the chiral limit. 

The whole approach can be extended to the case of non-degenerate quark masses, which is 
relevant for heavy flavour phenomenology. The only change is in the form of the WTI's. For 
the axial WTI, 2mq should be replaced by mi -|- m2, while for the vector identity, a term 
proportional to the mass difference should be added. 

The approach can also be extended to deal with more complicated composite operators, 
such as four-quark operators relevant to study the effective weak Hamiltonian, though the 
number of operators which must be included may render the procedure impractical. 

We stress again, however, that our approach fails, in general, for operators which can mix 
with lower dimensional ones (such as the octet part of the weak Hamiltonian). In this 
situation, lower dimensional operators may give rise in correlators to a large momentum 
behaviour that cannot be distinguished, and thus disentangled, from that coming from the 
exchange of Goldstone bosons in correlators with the insertion of the original operator. 

Appendix A 

In this Appendix we want to explain in a simple example why, in presence of spontaneous 
breaking of chiral symmetry, it may not be possible to fix non-perturbatively all the coeffi- 
cients of the mixing of lower dimensional operators with higher dimensional ones by simply 
looking at the high momentum behaviour of correlators. The reason is that correlators with 
the insertion of lower dimensional operators may turn out to have the same high momentum 
behaviour as that coming from contributions due to the exchange of Goldstone bosons to 
(some of the) correlators in which the original operator is inserted. 

To prove this statement let us consider the correlator [| 

GiM) ^ J d'xd'ye-'P''+'P'y{6A{q{y)q{x)O,{0))) , (68) 

where Od is an operator of dimension d and Sa{0) represents the axial variation of the 
operator O. A simple dimensional counting shows that the asymptotic scaling behaviour of 

^As in the rest of the paper also in this Appendix to simphfy notations we will neglect flavour indices. 
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G^{p,p') for large p and p' is given by the formula 



Grf(Ap,Ap') ~ A-^-^ (69) 
A similar analysis for the Green function 

G7ip,p') ^mj d'xdSd'ze-'P^+'P'y{q{z)^,q{z)q{y)q{x)OM) , (70) 
in which the operator P = q'j5q is inserted together with Od, gives 

G^{Xp,Xp') A'^-^ (71) 

A— +00 

Suppose now we are interested in the case of the mixing of an operator of dimension 6 (like 
the octet part of the effective weak Hamiltonian) with an operator of dimension 3 (like q'^^q 
or qq). From the scaling behaviours (|69D and ( ffTD we immediately conclude that, in the 
absence of a massless pion, everything works fine, because at large momenta 

1) the contribution of the correlator to the relevant axial WTI is subleading with respect 
to that coming from both for d = Q and d = 3; 

2) correlators with inserted operators of dimension 6 and 3 have well separated asymptotic 
behaviours. 

In this circumstance the unwanted contributions (i.e. those corresponding to (i = 3 in the 
above scaling laws) can be isolated in the various correlators and set to zero by appropriately 
choosing the mixing coefficients at our disposal. 

The problem arises if chirality is spontaneously broken, because in the chiral limit the pion 
contributes to GIT a term 



Gr(p,p') . = U d'xd^ye-^P^-^'P'y{n\q{y)q{x)Oem 

pion J 

J d'xd'ye-'P^+'P'ycwix,y){7T\q^,q{0)) , (72) 



where fjr is pion decay constant. The first equality in Eq. (^) follows from PCAC, the 
second from inserting the leading non-vanishing contribution of the O.P.E. 

q{y)q{x)OG{0) ~ Cw{x, y)qi^q{^) + ... . (73) 
A straightforward dimensional argument leads to the asymptotic scaling law 

cw{Xp,Xp')= I d^xdSe-'^P''+'^P'ycw{x,y) ~ A^^ . (74) 

The trouble with this behaviour is that it is the same as that predicted by Eq. ( |69D for 
an operator of dimension c? = 3. This simple observation proves our theorem. In fact, in 
the "kinematical" situation of the example we are considering, it would be impossible to 
distinguish at high momenta the contribution to the axial WTI coming from an operator of 
dimension 3 from the similar one due to the exchange of a massless pion in G^. 
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Appendix B 



In this appendix we wish to give a simple argument to show that in a gauge-fixed theory 
BRST symmetry does not forbid the mixing of gauge invariant operators with gauge non- 



invariant ones if the latter vanish by the equations of motion |^2|, P3 |. 

To simplify notations we collectively indicate by [0fc(a;)] the set of fundamental (fermionic 
and bosonic) fields, including ghosts, appearing in the Lagrangian of the theory and by 
their (bosonic and fermionic) BRST variations 0: 

5B0fc(a;) = Bk{x) . (75) 

Introducing the sources J^(a;), coupled to the variations Bk{x), we define the so-called Zinn- 
Justin functional by 

Z[j, Jb,v] = jVcP exp{-5[0] + j (j^kf + J J^BBk + J vO} = 

= jvct) exp{^[0, Jb] + / <Pk3^ + J vO} , (76) 

where r] is the source coupled to the BRST-invariant operator O and a summation over 
repeated indices is understood. Expectation values with respect to this measure will be 
denoted 

In the second equality of Eq. ( |75D we have introduced the definition 

A[(I),Jb] = -S[<P]+ J j'^Bk. (78) 

We assume that the operator O is renormalized, so that its insertions with the fundamental 
fields of the theory are finite if the fields and O are physically separated. The issue at hand 
is whether there can be contact terms, whose coefficients diverge when the regularization is 
removed, and the removal of which requires the use of gauge non-invariant operators. 

All the consequences of the BRST symmetry can be compactly expressed by the so-called 
Zinn- Justin equation 

/-^^■'<^'-- (-) 

This is derived by performing the change of variables induced by the BRST transforma- 
tions (^) in the functional integral ([76|) , assuming that the Jacobian is unity. It is useful to 
expose the condition on Bk that a unit Jacobian implies: 

= Jdxjv<f) — ^ f^Bkix) exp{^[0, Jb] + J + / 



dx 



6(f)k{x) 
SBk(x) 



+ ( / dxBk{x)jk{x)) . ^ (80) 



j,JB,V 



^Throughout this appendix to avoid cumbersome equations we will be rather cavalier about signs. 
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where we have used the BRST invariance of the action, 



= SbS = J dxBk{x 
and the nilpotency of BRST transformations, 



5S 



(81) 



= j dyB,{y) 



The Zinn- Justin equation (|79|) follows from 



6(f)k{x) 

as long as the condition for unit Jacobian 



821 



j,JB,V 



0, 



(83) 



is satisfied. One readily verifies that antisymmetry arguments imply that this equation is 
trivially satisfied. Indeed, the BRST variations, B^, are such that 



SBkjx) 
6(t)k{x) 







^4) 



We note that the Zinn- Justin equation holds in a regularized theory as long as regularized 
BRST transformations can be defined under which the regularized action is invariant, and 
which satisfy (|83|). This is the case for lattice QCD [^(see sec. 2.1). 

The regularized Green functions of the operator O obey the Slavnov- Taylor identities that 
follow from Eq. (^). If these Green functions contain contact terms which diverge when the 
regulator is removed, then these contact terms must separately satisfy the identities. Thus 
we can choose the counter-terms we add to cancel these divergences so that the identities, 
and thus the Zinn- Justin equation itself, are still satisfied. Thus we introduce the modified 
functional 



Z[3, Jb,v] = / exp{A[(P, Jb] + 



7^{0 + F[JbA])}. 



15) 



where F[Jb,(P] is a local functional of the sources, J^, and the fundamental fields, (pk, and 
derive the constraints imposed on the form of the functional F by the modified Zinn- Justin 
equation 



^ J2[j,Jb,v] k. . 



0. 



(86) 



The significance of the dependence of F on will become apparent in the concrete example 
we give below. 



^7) 



A little calculation shows that (|86D is equivalent to the equation 

5'F{y) 



dyr){y) 



dx 



5J^{x)5(pk{x] 



+ WF{y) 



+ 



+ / dyr]{y) / dy'riiy') / dx 



6F{y) 6{0{y') + F{y')) 
'SJ^{x) 6(f)k{x) 
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where we have introduced the operator 



W = dxi 1 

J ^6(j)k{x) 6J^{x) 6J^{x) 6(f)k{x] 



^6(j)k{x) 5J%{x) S(f)k{x) 

SO that 



To get Eq. (p7\} we have exploited the identity (|8^) and the BRST-invariance of the gener- 
ahzed action (Eq. ([75|)), which in our notation amounts to the two equations (^T]) and (|5^. 
Limiting for simphcity the analysis to Green functions with single insertions of the operator 
O, Eq. ( |87D implies the condition 

where we have reinstated the dependence on h, and defined the new operator W. 



The key observation needed to solve (^) is to recognize that is a nilpotent operator, 
W'^ = 0. Using this property, and the fact the fermionic operator / dy6'^ / 6J^{y)6(f)k{y) is 
nilpotent, it then follows that W itself is nilpotent, W'"^ = 0. Thus a solution to (0) is 
F = W'C. To show that this is the only solution, to all orders in perturbation theory, we 
expand F in powers oi%: F = Fq + hFi + h^F2 + h^F^ + . . .. Introducing this expansion in 
Eq. ( pOD and repeatedly using the nilpotency of W and W, one obtains the solution in the 
form^ 

F = wiCo + hCi + h^C2 + h^Cs + ...) + 

+^ / ^^^XT177^(^o + + ^'^2 + ...), (91) 
J 6(j)k{y)dJ^{y} 

where Ci[(f), Jb] are arbitrary local functionals of the sources, Jb, and the fundamental fields 
of the theory. The series expansion in h in the r.h.s. of Eq. (pl| ) can be formally resummed, 
giving 

F(.) = W'C(.) = WC(.) + H 1 1yj^§§^^ . (92) 

Remembering the form of W (see Eq. (^)), we can split F in a BRST-invariant and a non 
BRST-invariant part: 

F{x) = F^Bix) + 5BCix) (93) 
f 6A 6C(x) f S'^Cix) 



^To begin the iteration one needs the result that solutions to WG — are of the form G — WH, where 
iJ is a local functional of fields and sources fS^l . 



20 



We are really interested in -F/vb, because the BRST-invariant part of F can always be 
reabsorbed in a redefinition of O. We can now see that, as expected, the insertion of F/vb 
with the fundamental fields of the theory gives rise only to contact terms, thanks to the 
locality of the functional C. To see this we consider the Green function 



Inserting Eq. 



one gets, after an integration by parts, 
G{x, xi, 



(95) 



(96) 



i=l 



where we have used the definition 

6C{x) 



SJB'iXi) 



XkXx)S{x - Xi 



(97) 



If we now set = in Eq. ( |96D we see explicitly that -F/vs generates contact terms, and 
that they are not constrained by BRST symmetry (since C and thus the are arbitrary). 

Returning to the full expression of F, we can now rewrite Eq. (p2|) as 



F{x) 



-Xfc(x) + SbC{x) +h- ^ 



(9^ 



Apart from the last term, Eq. (p8| ) proves our initial statement, i.e. F vanishes by the 
equations of motion up to BRST-invariant terms. The last term is proportional to 5(0). It 
is present to ensure that the term vanishing by the equations of motion does in fact lead to 
a contact term [as seen in the derivation of Eq. (p6D above] . More formally, it compensates 
divergent contact terms (hidden) in the product of operators appearing in the first term of 
Eq. (pHl). For the flavor non-singlet operators considered in this work, one can in fact show 
that this term is algebraically zero. 

We want to end this appendix by clarifying, in a significant example, the physical meaning 
of the presence of the counter-term F in Eq. (^). For this purpose let us consider the 
Slavnov- Taylor identity 

+ 



5j"(z)5J]^(t)5r/(a;) 6j"'{t)6J'^{z)6T]{x) 







(99) 



j=Jg=rj=0 



that follows from equation (|86|) by successively differentiating with respect to jm(t), jn{z) 
and rj{x). After some algebra one gets 

^°To avoid introducing a new symbol we now use Xk to indicate the functional derivative 6C /SJ^ evaluated 
at Jb — 0. 
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+6{X - Z) (Xn{z)Bra{t)) + - t) (0n (^)<5b (X^ (t) ) ) + 

- Z)b{x - t) {Xnm{x)) 



sMi.m.) ^„ - iMm) + (101) 

+5{x - t) {X^{t)B^{z)) + 6{x - z) {(j)^{t)6B{Xn{z))) + 

+6{X - Z)6{X - t) {Xmn{x)) . 

In the above equations we have used the formal expansion |^ 

oo 

C{x) = Y: Jb'{x) . . . j'dix)X,,...k„{x) . (102) 

n=l 

In Eqs. ( |100| ) and ( |101| ) we see explicitly the meaning of the terms in Fnb proportional to 
Jb (which arise in part from terms in C quadratic in Jb)'- they lead to the double contact 
terms proportional to Xnm- Notice that the sum of all contact terms adds up to zero by 
itself by virtue of the (anti) symmetry properties of Xnm and of the BRST identity obeyed 
by {Bn{z)Xm(t)) , which reads 

{XUt)Bn{z)) = {5B{Mz))XUt)) = - {M^)SB{XUt))) ■ (103) 

The cancellation of contact terms is also consistent with the equation 

{Bm{t)Mz)0{x)) + {<f)m{t)Bn{z)0{x)) = {6B{Mt)M^)0{x))) = (104) 

that follows from the BRST invariance of O. 

Let us summarize what we have learned. It is consistent with BRST invariance for the gauge- 
invariant operator to have (divergent) contact terms with the fundamental fields. These can 
be cancelled by appropriate choices of the X„ and Xnm in such a way that the resulting 
correlation functions have finite Fourier transforms. In other words, adding the counter- 
term F amounts to redefining the correlation functions as follows: 

{<j)m{t)Bn{z)0{x)y = {<Pm{t)Bn{z)0{x)) + 

+5{X - Z) {(l)m{t)5B{Xn{z))) + 5{x - t) (X„(t)E„(t)) + 

+5{x-z)5{x-t){Xmn{x)) , (105) 
SO that they are finite for all x, z and t. The Slavnov- Taylor identity (pOf ) then implies that 

{cPm{t)Bn{z)0{x))' + {cPn{z)Bm{t)0{x))' = , (106) 

i.e. {(pmit) Bn{z)0 {x))' and {(j)m{t) Bn{z)0 (x)) satisfy the same BRST identity. 

^^The term with n = is dropped from the expansion because a Jg-independent piece would give a 
BRST-invariant contribution to O. 
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What are the constraints on the counter-term imphed by BRST invariance? Although X„, 
the contact term between 0„ and O, is not constrained, the contact term between and 
O is required to be This relationship between contact terms is encoded in the Jb 

dependence of F. If more than two fields collide, then additional arbitrary contact terms 
are allowed. This way of understanding the result makes it clear that the general form of 
the contact terms could, with a lot of imagination, have been guessed without reference to 
the Zinn- Justin equation. 



Appendix C 

In this appendix we test our method for off-shell improvement using one-loop perturbation 
theory. We focus on the quark propagator, the improvement of which requires, we claim, 
the gauge non- invariant term proportional to cngi- The necessary one- loop perturbative 
calculations have been worked out by Capitani et al, and presented in final form in Ref . ||2^ . 



We stress, however, that the method of off-shell improvement presented in Ref. ||2^ does not 
include the c^gi term. It works, however, up to one-loop, because, as we shall show in this 
Appendix, cngi actually vanishes at tree-level. 

It is convenient to work with the inverse bare lattice propagator since this is the quantity 
calculated in perturbation theory. Thus we rewrite Eq. ( P?] ) as 

= Sz\p) (l + ab,m - 2ad,S-L\p) - 2ac^Gii /+ 0{a^)) . (107) 

Here we have used the fact that [S£^,;^ = (since the only gamma-matrix structure allowed 
in Sj^^ is ]/), and expanded the quark- field normalization factor in the conventional way 

Zg = Zj(l + abgm) . (108) 

The advantage of moving Z^ to the l.h.s. of Eq. (|107|) is that it separates the issues of 
improvement and normalization. The l.h.s. is proportional, up to 0{a^), to the continuum 
inverse propagator, and thus it should be possible to choose bg, c'g and c^gi so that the r.h.s. 
does not contain 0{a) terms. An important constraint is provided by the fact that these 
three improvement coefficients should be functions only of the bare coupling (7o-0 

The one-loop result for the inverse propagator with an improved action takes the form 

2 

Sl\p) = 2]/+m + ^-ASi + 0(a2), (109) 
El = i + mi:2L + ap'^^sL + ai ^mll^L + arn^Tj^L -\ 5— SeL , (HO) 
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Strictly speaking, improvement coefRcients are functions of the coupling g^ll + bgam), as explained in 
Ref. j7|. This distinction is, however, important only at two-loop order for the quark propagator. 
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where explicit results for the S^x (for arbitrary m? /p^) can be deduced from the results in 
Ref. [^. Here the quark mass is defined as usual 



using the one-loop result for Kc. We need only the large p behavior of the S^x, 

Eji = £j In(aV) + dj + O(mVp^) • (113) 

The coefficients of the logarithm are (using the appropriate leading order value, csw = 1) 

4 =£3 = a, £2 = -^4 = -24 = (3 + «), 4 = 0, (114) 

where a is the gauge parameter which vanishes in Landau gauge. It is important to note 
that subleading terms proportional to {m? /p^) are absent from Ei^ and Tj^l-, since these are 
explicitly included in SgL and respectively. These terms do contribute in the large 
limit. However, the other such terms (in S2L and — Sg^) give vanishing contributions 
in this limit .0 

Using these results, we can evaluate the r.h.s. of Eq. ( p.07| ) 



= z y(l - X^il) + m(l - AS2l) 
+ a/ [2c;(l - 2ASil) + 2cNGi{l - ASu.) + 1/2 - ASsl] 

+ ai ]/m [-44(1 - ASii - XL^l) - 2cngi{1 - ASsl) + K{1 - ASi^) - XL^l] (115) 
+ am^ [-24(1 - 2AS2L + 2AS6l) + 6,(1 - AS2l) - cngi'2XT,^l - ASjl] + 0(a^ A') . 

From this we can deduce the constraints on the improvement coefficients. First we note that, 
for large momenta, 

S'\p) = S-\p) + 0{^) = --l^^ fpy) + Oiml/p-) , (116) 

where we have used Eqs. ( p3|) and ([26|) . Asymptotically, Si(j9^) is independent of rrig, while 
S2/S1 = niq up to logarithmic corrections in p^. It follows that, asymptotically, there are no 
terms proportional to p^ and i jkn in 5*"^, and so the coefficients of these terms in Eq. (|115|) 
must vanish. These constraints can be written (up to corrections of 0(A^)) 

24 + 2CNGI + 1/2 = A [(44 + 2CNGI + l)SiL + (S3L - , (117) 

44 + 2c^G/ - 6g = A [(44 - 6,)Sii + (44 + 2cjvG/ + 1)S2l - (S2L + ^4^)] • (118) 

The final constraint is obtained by noting that S^^ depends linearly on the renormalized 
mass rriq, which is related to the lattice quark mass by 

rriq = ZmTn{l + abmin) . (119) 



""^•^We thank Paul Rakow for pointing out to us the importance of including the apparently subleading form 
factor SgL, as has been done in Ref. | |27[ |. 
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It follows that the ratio of and m terms in Eq. ( p.l5| ) gives bm- In this way one finds 

+ 2c; - 6, = A [(2c; + 1/2)S2L - (4c; + 2cNGi)^6L - (S2l/2 + Ssl)] + 0{\^) . (120) 

The three equations ([TT7|) , ([TTsD and (|l2g) provide non-trivial constraints because their 
l.h.s.'s are constants, independent of ap, while the r.h.s.'s contain terms with logarithmic 
dependence on ap. At first sight, the system is over-constrained, because each of the three 
equations has an 0(1), an 0{gQ), and 0{gQ\n{ap)) part, providing nine constraints in total, 
while there are four improvement constants, each with an 0(1) and an 0((?o) part to be 
determined. It turns out, however, that the constraints are redundant, and, in fact, the 
system is under-constrained. 

To proceed we expand the improvement constants in powers of A, e.g.[3 

6^ = 6(o) + A6«+0(A2). (121) 

Then the first constraints are that the l.h.s.'s of the three equations vanish at leading order, 
since the r.h.s.'s are of 0(A). This gives 

2c;(°) + 2cSL/ + 1/2 = 0, 4cf^ + 2cSL/ - = , b^^^ + 2c',^'^ - bf = . (122) 

We see that, as noted in the text, we cannot determine the tree- level value of all four constants 
with a tree- level computation of the propagator. Only the combinations c'q^^'^ + c^^^qj = —1/4, 
2c; - = 1/2 and b^^^ = -1/2 are determined. 

The one-loop computation does, however, completely determine the tree-level improvement 
coefficients. To see this, note that we have written the three equations so that the combina- 
tion of in the last term on each r.h.s. has no dependence on In(a^p^) [see Eq. ( |114| )], and 
is thus just a constant (in any covariant gauge). In Landau gauge, Si^ is also a constant, 
but has a logarithmic part. To cancel this logarithmic part on the r.h.s.'s of Eqs. (|118|) 
and ( |120| ) requires 

4c; + 2cS^, + 1 = 0, 2c;(°) + 1/2 = 0. (123) 

These two extra equations could be inconsistent with the previous three, and is a non- 
trivial test of improvement that they are not. They provide the extra constraint needed to 
determine the tree level improvement coefficients, and one finds that 

cPgi = . (124) 

Thus the "new" constant does vanish at tree level, but from the point of view of the cal- 
culation this is an accident. For the tree- level values of the other two constants one finds 
c;(o) = -1/4 and 6f = 1. 

We note in passing that the coefficient of T,il in Eq. (|1 1 7|) does vanish due to Eq. (|123|) . 
This means that the improvement works in perturbation theory at one-loop in any covariant 
gauge. 

^''Note that it is standard to expand in powers of rather than A. We choose to expand in A so as to 
simphfy equations. 
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Finally, we can deduce the constraints from the non-logarithmic 0(A) terms. These can be 
written 

2d,^^^ + 2c%i = S3L-Sii, (125) 

2c;(')-6« = -SsL + - - , (126) 

6« = SsL - + S2L/2 + - SsL + . (127) 

Thus a one-loop computation of the propagator does not determine c^-ngi- To do so would 
require a two-loop computation of the propagator. Alternatively, noting that the gauge non- 
invariant term can be rewritten as ^g, it seems likely that a one-loop computation of the 
quark- gluon vertex would determine c^}^gi- 
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